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ABSTRACT 

We show that recent measurements of the power spectrum of cosmic microwave background 
anisotropics by BOOMERanG and MAXIMA can be mainly characterized by four observables, the 
position of the first acoustic peak l\ = 206 ± 6, the height of the first peak relative to COBE normal- 
ization Hi = 7.6 ± 1.4, the height of the second peak relative to the first H2 = 0.38 ± 0.04, and the 
height of the third peak relative to the first H3 = 0.43 ± 0.07. This phenomenological representation of 
the measurements complements more detailed likelihood analyses in multidimensional parameter space, 
clarifying the dependence on prior assumptions and the specific aspects of the data leading to the con- 
straints. We illustrate their use in the flat ACDM family of models where we find £l m h 3 8 > 0.079 (or 
nearly equivalently, the age of the universe to < 13— 14Gyr) from l\ and a baryon density Vl^h 2 > 0.019, 
a matter density O m /i 2 < 0.42 and tilt n > 0.85 from the peak heights (95% CL). With the aid of several 
external constraints, notably nucleosynthesis, the age of the universe and the cluster abundance and 
baryon fraction, we construct the allowed region in the (Q, m ,h) plane; it points to high h (0.6 < h < 0.9) 
and moderate VL m (0.25 < Q, m < 0.6). 

Subject headings: cosmic microwave background — cosmology: theory 



1. INTRODUCTION 



Wi th the data from the BOOMERanG ( |de Bcrnardis 
ct al.| 2000) and MAXIMA ( [Hanany et"al| 2000) experi 
mcnts, the promise of measuring cosmological parameters 
from the power spectrum of anisotropics in the cosmic mi- 
crowave background (CMB) has come substantially closer 
to being fulfilled. Together they determine the location of 
the first peak precisely and constrain the amplitude of the 
p owen at the expected position of the accond peak. — The 



MAX IMA experiment also limits the power around the 



expected rise to the third peak. 

These observations strongly constrain cosmological pa- 
rameters as has been shown through likelihood analyses 



prior assumptions (|Lange et al. 


2000; |Balbi et al. 


Tegmark & Zaldarriaga 2000b; 


Bridle et al. 2000). 



While 



dimension of the parameter space makes it difficult to un- 
derstand what characteristics of the observations or prior 
assumptions are driving the constraints. For instance, it 
has been claimed that the BOOMERanG data favor closed 



universes ( White ct al. 20 00; Langc et al. 2000) a nd high 
baryon density (|Hu[2000; [Tegmark fc Zaldarriaga] 2000b), 
but the role of priors, notably from the Hubble constant 
and big-bang nucleosynthesis is less clear. Indeed, whether 
CMB constraints agree with those from other cosmological 
observations serves as a fundamental test of the underly- 
ing adiabatic cold dark matter (CDM) model of structure 
formation. 

In this paper we show that most of the information in 
the power spectrum from these two data sets can be com- 
pressed into four observables. The correlation among cos- 
mological parameters can be understood by studying their 
effects on the four observables. They can also be used to 
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search for solu tions outside the sta ndard model space (e.g. 
Peebles et al.| 2000; |Bouchet et alj 2000). 



AT an instructive application of this approach, we con- 
sider the space of flat adiabatic CDM models. Approx 
imate flatness is clearly favored by both BOOMERanG 
and MAXIMA (|de Bernardis et alj 2000; |Hanany et al 



2000) as we ll as previou s data, notably from the TOCO 



exper iment ( Miller ct al. 1999), as shown by pr evious anal- 
yses ( | Lincwoavci 1998; pfstathiou ct al. 1999; Tegmark & 



Zaldarriagaj 2000a). 

Our main objective in this application is to clarify the 
constraints derived from the CMB observations using the 
likelihood analyses and understand how they might change 
as the data evolves. Then, with the aid of a few external 
constraints we map out the allowed region in the plane of 
the matter density (f2 m ) versus the Hubble constant (Ho: 
we use h to denote the Hubble constant Ho = lOO/i km 



L Mpc~ 



The external constraints which we employ 



include (i) the rich cluster abundance at z « 0, (ii) the 
cluster baryon fraction, (iii) the baryon abundance from 
big bang nucleosynthesis (BBN), and (iv) the minimum 
age of the universe. We also discuss their consistency with 
other constraints, such as direct determinations of Hq, f2 m , 
and the luminosity distance to high redshift supernovae. 
All errors we quote in this paper are at 67% confidence, 
but we consider all constraints at a 95% confidence level. 

In §2 we start with a statistical analysis of the CMB 
data. We introduce the four observables and discuss their 
cosmological implications. In §3, we place constraints on 
the (Q m ,h) plane and discuss consistency checks. In §4, 
we identify opportunities for future consistency checks and 
arenas for future confrontations with data. We conclude 
in §5. The appendix presents convenient formulae that 
quantify the cosmological parameter dependence of our 
four characteristic observables in adiabatic CDM models. 
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Fig. 1. — Power spectrum data and models: (left panel) full range on a log scale; (right panel) first 3 peaks on a linear scale. The 
BOOMERanG (BOOM) and MAXIMA (MAX) points have been shifted by their lcr calibration errors, 10% up and 4% down respectively. 
Dashed lines represent a model that is a good fit to the CMB data alone: Q m = 0.3, Oa = 0.7, h = 0.9, fij/i 2 = 0.03, n = 1 which gives 
l\ = 205, Hi = 6.6, H2 = 0.37, H3 = 0.52. Solid lines represent a model that is allowed by our joint constraints: Q m = 0.35, Q\ = 0.65, 
h = 0.75, fl b h 2 = 0.023, n = 0.95 which gives £1 = 209, H 1 = 5.8, H 2 = 0.45, H 3 = 0.5. Note that the labelling of the H's in the figure is 
schematic; these values are the power ratios as defined in the text. 



2. CMB OBSERVABLES 

2.1. Statistical Tests 

With the present precision of the BOOMERanG and 
MAXIMA observations (see Fig. [I]), it is appropriate to 
characterize the power spectrum with four numbers: the 
position of the first peak £1, the height of the first peak 
relative to the power at £ — 10 



VAT10 



the height of the second peak relative to the first 



H 9 



A7> 



and the height of the third peak relative to the first 



H* = 



AT £; 
AT,. 



(1) 



(2) 



(3) 



where (ATi) 2 = £(£+l)Ct/2it with C, the power spectrum 
of the multipole moments of the temperature field. Note 
that the locations of the second and third peaks are set 
by t heir harmonic relation to the first peak [see Appendix, 
eq. (A7)] and so H 2 and H 3 are well-defined even in the 
absence of clear detections of the secondary peaks. 

One could imagine two different approaches towards 
measuring these four numbers. We could extract them 
using some form of parametrized fit such as a parabolic fit 



to the data ( |Knox fc Page] 2000; de Bernardis et al.| 2000). 
Alternately, we could u se template CDM mo dels as cal- 



culated by CMBFAST ([Beljak fc Zaldarriagaj 1996), and 
label them by the values of the four observables. We can 
measure \ 2 f° r these CDM models and interpret them as 
constraints in the four observables. Both of these methods 



give similar results. We chose the second one because it 
is more stable to changes in the £ ranges taken to corre- 
spond to each peak and incorporates the correct shape of 
the power spectra for CDM-like models. 

To determine the position of the first peak, we take the 
data that fall between 75 < £ < 375 and carry out a \ 2 
fitting using a flat model template (O ro + = 1 here 
and below unless otherwise stated) with varying h and 
Vt m at the fixed baryon density flbh 2 = 0.02 and tilt pa- 
rameter n = 1. We include calibration errors, 10% for 
BOOMERanG and 4% for MAXIMA. Figure | shows A X 2 
as a function of £\ for the BOOMERanG data alone and 
for the combination of BOOMERanG and MAXIMA. The 
figure implies 



£i= 200±8 (BOOM), 

£1 = 206 ± 6 (BOOM + MAX) . 



(4) 



Other choices of Sib ft 2 and n for the template parameters 
modify slightly the value of \ 2 but not A\ 2 or the al- 
lowed region for l\. It is noteworthy that adding in the 
MAXIMA data steepens Ax 2 on both sides of the mini- 
mum despit e the preferenc e for £\ ~ 220 in the MAXIMA 
data alone (Hanany et al. 2000). The fact that both data 
sets consistently indicate a sharp fall in power at I > 220 
increases the confidence level at which a high £\ can be 
rejected. 

The Hi statistic depends both on the acoustic physics 
that determines the first peak and other processes relevant 
at £ ~ 10. The shape of the template is therefore more 
susceptible to model parameters. We choose to vary n 
which changes Hi and also allow variations in £!a so that 
the position of the first peak can be properly adjusted. 
The other parameters were chosen to be Vl^ti 1 = 0.03 and 
n m h 2 = 0.2. Using the BOOMERanG and MAXIMA 
data for 75 < £ < 375 in conjunction with the COBE 
data, we find 

#1=7.6 ±1.4. (5) 



3 




A BOOM 

a BOOM+MAX 



200 



250 



h 



300 



350 



Fig. 2. — Constraints on the first peak position: Ax 2 (^i) for the 
data from 75 < £ < 375. We define the Ict errors to be 1/2.5 of the 
errors at 2.5ff (Ax 2 = 6.2, solid lines). 
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Fig. 4. — Constraints on the height of the third peak relative to the 
first: Ax 2 (H 3 ) for the data fro^ 75 < £ < 375 and 600 < t < 800. 
lcr errors are defined as in Fig. 



Other template choices can modify the constraint slightly 
but the errors are dominated b y the COBE 7% cosmic 
variance errors ( |Btmn k White| 1997) and MAXIMA 4% 
calibration errors on the temperature fluctuations. 
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Fig. 3. — Constraint on the height of the second peak relative to 
the first: Ax 2 (Hal for the data from 75 < £ < 600. Iff errors are 
defined as in Fig 

For H 2 , we take the data for 75 < t < 600 and con- 
sider templates from models with varying flbh 2 , n and 
r^A, where the last parameter is included to ensure that 
the models reproduce the position of the first peak. In this 
case, x 2 as a function of Hi minimized over exhibits 
some scatter due to information that is not contained in 
the ratio of the peak heights (see Fig. 0). Nonetheless, 
the steep dependence of A% 2 on Hi indicates that this 
statistic is robustly constrained against the variation of 
the template. Taking the outer envelope of A% 2 , we ob- 
tain 

H 2 = 0.37 ±0.04 (BOOM), 

Hi = 0.38 ± 0.04 (BOOM + MAX). (6) 

7?3 is only weakly constrained by the two highest I 
points (600 < t < 800) from MAXIMA in conjunction 



with the first peak data (75 < I < 375) from both exper- 
iments. We consider templates from models with varying 
Q m h , n, £lf,h 2 and The latter two parameters are 
included to ensure that the position of the first peak and 
the depth of the first trough can be modeled. Minimizing 
X 2 over these two parameters, we plot Ax 2 as a function 
of H 3 (see Fig. |) to obtain the bound 



H 3 = 0.43 ± 0.07 (BOOM + MAX). 



(7) 



Note that the constraints on H 3 employ a template-based 
extrapolation: points on the rise to the third peak are used 
to infer its height. 
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Fig. 5. — Peak position constraint in the (f2 m , C2a) plane with 
various priors. The prior assumptions weaken from the light shaded 
region to the dark shaded region and means that the unlisted 
priors are unchanged from the neighboring region of stronger priors. 
The constrained region is strongly limited by the range in h consid- 
ered and to a lesser extent that in Q b h 2 and the equation of state 
of A, w. 
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2.2. Cosmological Implications 

The values for the four observables we obtained above 
can be used to derive and understand constraints on cos- 
mological parameters from the experiments. 

The position of the first peak as measured by l\ is de- 
termined by the ratio of the comoving angular diameter 
distance to the las t scattering epo ch and the sound hori- 



zon at that epoch (IIu fc Sugiyamo 1005). Therefore it is 



para meter that depends only o n ge ometry and sound wave 

dynamics [see Appendix, eq. ( [A3] )] through f2 TO + h, 
i7 TO , u>, and Q, b h 2 , in decreasing order of importance. The 
equation of state parameter w = pa/ PA with pa and pa 
the pressure and energy density of the vacuum (w = —1) 
or negative-pressure energy; we use A to re fer to either 
option. The effect of tilt is small [see eq. (f 



Ali 
IT 



0.17(n- 1) 



(8) 



and so we neglect it when considering models with n ~ 1. 

Figure || displays the constraints in the O m — J7a plane. 
Notice that the confidence region is determined not by 
uncertainties in the measured value of £ 1; but rather the 
prior a ssumptions a bout t he acceptable range in fo, fl b h 2 
and w ( |Lange et al. 2000; Tegmark & Zaldarriaga 2000b). 
Given the broad consistency of the data with flat models 
(Oa + = 1) with w — — 1, we will hereafter restrict 
ourselves to this class of models unless otherwise stated. 




Fig. 6. — Peak positions in the (fi m ,/i) plane for Cli,h 2 = 0.019 
and 0.025. Allowed region is indicated by hatching and assumes 
tt h h 2 > 0.019. 



To better understand the dependence on the Hubble 
constant, we plot in Figure |] contours of constant t\ in 
the (Q. m ,h) plane for Vl b h 2 = 0.019 and Q b h 2 = 0.025. A 
higher baryon abundance decreases the sound horizon at 
last scattering and pushes up the contours in the direction 
of higher f2 m , h. 

The 95% limit, t\ < 218 from equation (||), excludes 
the lower left region shown in Figure H. This limit derived 



from the curve for il b h 2 = 0.019 is robust in the sense that 
this baryon abundance represents the mimimum value al- 
lowed by the CMB, as we shall see later. This constraint 
is summarized approximately by £l m h 3 S > 0.079. The 
dependence differs from the familiar combination fl m h 2 
since in a flat universe part of the effect of lowering f2 m is 
compensated for by the raising of Q\. The lower limit on 
fl m h 2 or equivalently on fl m found in Tegmark & Zaldar 



riaga (2000b) reflects the prior upper limit on h, e.. 
h < 0.8, n m h 2 > 0.12 or 0, 



for 



> 0.18. 

The 95% confidence lower limit of the combined fit is 
l\ > 194, but it does not give a robust limit in this plane 
in the absence of an upper limit on fl b h 2 . The lower limit is 
also less robust in the sense that it comes mainly from the 
MAXIMA result whereas both experiments agree on the 
upper limit in the sense that the addition of MAXIMA 
only serves to enhance the confidence at which we may 
regect larger values. 




200 



220 



h 



240 



260 



Fig. 7. — Age-^i correlation. Points represent models in the 
(£lm,h) plane. Jjines represent the simple fit to the outer envelope 
given in eqn. (0) for the two choices of f2;,/i 2 . Light shaded region 
represents current constraints on i\ ; dark shaded region is also con- 
sistent with to > 11 Gyr. 

For a flat geometry, the position of the first peak is 
strongly correlated with the age of the universe. The cor- 
relation is accidental since i\ is the ratio of the conformal 
ages J dt/a at last scattering and the present that enters 
not simply the physical age today. Nonetheless, Fig 
shows that the correlation is tight across the (f2 m , h) va' 
ues of interest. The upper envelope corresponds to the 
lowest Q m (= 0.1) and implies 



lGyr 



< 19.7 -0.155(250 -4) (0.68) 



l+u) 



n b h 2 

0.019 



1.7 



(9) 



where we have included a weak scaling with w. The obser- 
vations imply < < 13- 14 Gyr if fl b h 2 > 0.019. While this 
is a weak constraint given the current observational uncer- 
tainties, notice that the central value of the BOOMERanG 
results t\ ~ 200 would imply to « 9 — llGyrs. 

The H2 statistic, the ratio of the heights of the second 
peak to the first, mainly depends on the tilt parameter and 
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Fig. 8. — Peak height ratio H2 in the (n,f2(,/i 2 ) plane for 
Q m h 2 = 0.15 and 0.20. Allowed region is indicated by light shad- 
ing and assumes C m /i 2 > 0.15. The constraint at n ss 1 is nearly 
independent of this assumption. Dark shading indicates the region 
consistent with Hi (n > 0.85) and nucleosynthesis (fTl^h 2 < 0.028). 



the baryon abundance. This combination is insensitive to 
reionization, the presence of tensor modes or any effects 
that are confined to the lowest multipolcs. The remaining 
sensitivity is to fl m h 2 and is modest in a flat unive rse du e 
to the cancellation of two effects [see Appendix eq. ( A19 )]. 
Figure || shows contours of H2 =const. in the (n, Vt^h 2 ) 
plane for Sl m h 2 = 0.15 and 0.20. 

The result from the previous subsection, H2 < 0.46 at 
a 95% confidence level is shown with shades, giving a con- 
straint 



n h h 2 > 0.029 



tt h 2 
0.15 



-0.58 



(n- 1) 4-0.024. (10) 



Note that the constraint at n — 1 is approximately in- 
dependent of Vl m h 2 . This limit agrees well with the those 
from the detailed likelihood analysis of Fegmark & Zaldar- 
riaga[](2000b, c.f. their Fig. 4) under the same assump- 
tions for fl m h 2 and supports the claim that H 2 captures 



most pf the information from the data on these paramc- 
/Vn upper limit from H2 > 0.32 also exists but is 



ters. 

weaker than conservative constraints from nucleosynthesis 
as we discuss in the next section. 

We can derive a limit on n from the indicator H±. The 
cosmological parameter dependence of Hi is more compli- 
cated than the other two we discussed above. Fortunately, 
most complications tend to decrease Hi by adding large 
angle anisotropics. The lower limit on n from the lower 
limit on Hi in the absence of e.g. reionization or tensor 
modes is therefore conservative. The upper limit on n 
is very weak unless one excludes the possibil ity of tensor 
modes as a prior assumption (c.f. Balbi et al. ). We search 
for the minimum n that gives Hi larger than the 2a lower 
limit Hi > 4.8 along the parameter space that maximizes 
Hi- This gives a conservative lower bound of n > 0.85. 
This bound is to good approximation independent of £l m h 2 
for Sl m h 2 > 0.15. Below this value, the bound tightens 



marginally due to the integrated Sachs- Wolfe (ISW) effect 
on COBE scales, but in such a way so as to maintain the 
bound ilbh 2 > 0.019 when combined with th e constraint 
from H2, th e inequality (Eoh. The analysis of Tcgmark & 



Zaldarriaga (2000b) yields n > 0.87 in the same parame- 



ter and data space indicating that not much information 
is lost in our much cruder parameterization. 




Fig. 9. — Peak height ratio H3 in the (n,il m h 2 ) plane for Hi = 
0.38. Allowed region is indicated by light shading. Darker shading 
indicates the allowed region with additional constraints from Hi 
(n > 0.85) and nucleosynthesis (f2j,/i 2 < 0.028, which gives n < 
1.16). 

The -H3 statistic depends more strongly on fl m h 2 and 
n since the baryons affect the height of the third and first 
peak similarly. In Fig. ^, we show the constraint in the 
(n,fl m h 2 ) plane with the baryon density fixed by requiring 
H2 = 0.38. When combined with the constraint on the tilt 
n > 0.85, we obtain Q, m h 2 < 0.42. 

In Fig. |l| (dashed lines), we compare a model designed 
to have acceptable values for £%, Hi, H2 and H3 with 
the power spectrum data from COBE, BOOMERanG and 
MAXIMA. This gives \ 2 = 27.2 for 30 data points and 



compare s well with the best fit model of Tegmark & Zal- 



darriaga (2000b) with their "inflation prior" where \ z 



26.7. We summarize the constraints from the CMB as: 
n m h 3 8 > 0.079 (or t < 14 Gyr); n > 0.85; n b h 2 > 0.019; 
n m h 2 < 0.42. 



3. EXTERNAL CONSTRAINTS 

In this section, we combine the constraints from the 
CMB with those from four other observations: the light 
element abundances as interpreted by BBN theory, the 
present-day cluster abundance, the cluster baryon frac- 
tion and age of the universe. We then translate these 
constraints onto the (Q m ,h) plane and discuss consistency 
checks. 

3.1. Nucleosynthesis 

The first external constraint wc consider is that on 
baryon abundance from primordial nucleosynthesis (see 



6 



Copi et al. 1995; Olive et al. 1999; Tytler et al. 2000 



for recent reviews). |01ivc et al. (1999) give a high baryon 
density option of 0.015 < D.bh' 2 < 0.023, and a low baryon 
density option 0.004 < £l b h 2 < 0.010 as a 2a range. A 
low baryon density is indicated by the traditiona l low 
value of helium abundance (Y p — 0.234 ± 0.003) ( plivc 
et alj]l999), and agrees with a literal interpretation of the 
lithium abundance. There are also two Lyman limit sys- 
tems which taken at face v alue p oint to high deu terium 



abundance (Songaila et aj| 1994; Carswell et al. 1994; 



Buries et al 



1999a; Tytler et al 



1999) and implies a 

low baryon density. 

Our lower limit on £l b h 2 > 0.019 from Hi and H\ is 
strongly inconsistent with the low baryon abundance op- 
tion. In fact our limit is only marginally consistent with 
even the high baryon option if we take the latest deter- 
mination of the deuterium abundance at face value and 
treat the individual errors on the systems as statistical: 
D/H= (3.4 ± .25) x 10~ 5 from 3 Lyman limit systems 
dKirkman et al.| 2000 ) which impli es n b h 2 = 0.019±0.0012 



(Tytler et al, 2000; |Burles et al. 199 9b). In this paper 



we provisionally accept the 2a limit of Olive et al. (1999). 
The issue, however, is clearly a matter of systematic errors, 
and we discuss in what follows where they could appear, 
trying to find a conservative upper limit. 

Since it is possible that the measured D/H abundance 
is high due to contamination by H, we consider the firm 
lower limit on the D/H abundance fr om interstellar clouds. 
The earlier UV data (McCullough 1992) show a varia- 
tion of D/H from 1.2 to 2.5 x 10~ 5 . This variation is 
confirmed by modern high resolution spectrographs. The 
clouds studied are still few in numb er and range from 
D/H=(1.5±0.1) x 10~ 5 flLmskyl 1998; [Linsky et at] 1995) 
to 0.7 x 10 -5 ( penkins et al.| 1999). This variation is 
reasonable since the clouds are contaminated by heavy el- 
ements, indicating significant astration effects. Therefore, 
we take the upper value as the observational D/H abun- 
dance, and take the minimu m ast ration effect (facto r 1.5) 



from model calculations (see Tosi 1996; Olive et al. 1999, 



modified for a 10 Gyr disk age) to infer the lower limit on 
the primordial deuterium abundance. We take 2 x 10 -5 as 
a conservative lower limit on D/H. This value agrees with 
the p re-solar system de uterium abundance inferred from 
3 He ( ploeckler fc Geiss| 1998). This D/H corresponds to 
Vt b h 2 = 0.028 and Y p = 0.250 - 0.252. 

The helium abundance Y p directly depends on the the- 
oretical calculation of the helium recombination line, and 
the discrepancy between the estimate (Y p = 0.244 ±0.002, 



Izotov & Thuar 1998) and the traditional estimates (Y p 



0.234 ± 0.003) largely arise s from the tw o different calcu- 
lations, |Smits| (1996) and |Blocklehurst] (1972). The he- 
lium abundances derived from three recombination lines 
He I A4471, A5876 and A6678 for a given HII region differ 
fractionally by a few percent. Also, the effect of under- 



lying stellar absorption by hot stars is unclear: Izotov & 
Thuai§ (1998) use the departure of the He I A6678/A5876 
strengths from the Smits calculation as an estimator, but 
a calculation is not actually available for the He line ab- 
sorption effect. While these variations are usually included 
as random errors in the nucleosynthesis literature, we sus- 
pect that the error in the helium abundance is dominated 
by systematics and a further change by a few percent in 
excess of the quoted range is not excluded. 



The interpretation of the Li abundance rests on a sim- 
plistic model of stars. It seems that our understanding 
of the 7 Li abundance evolution is still far from complete: 
for instance we do not understand the temperature gradi- 
ent of the Li/H ratio in halo dwarfs, which shows a trend 
opposite to what is expected with 7 Li destruction due to 
diffusion. Hence we do not view the primordial 7 Li abun- 
dance determinations as rock-solid. 

We therefore consider two cases Vt b h 2 < 0.023 as a 
widely accepted upper limit and fl b h 2 < 0.028 as a very 
conservative upper limit based on interstellar deuterium. 
When combined with the limit of eq. ([l0|), the latter con- 
straint becomes n < 1.16 for £l m h 2 < 0.2 (as appropriate 
for setting a lower bound on f2 m in the next section, see 
also Fig. 0); if we instead take Q b h 2 < 0.023 (Olive et al 



1999), the limit becomes n < 0.98. In conjunction with the 
constraint from Hi , the allowed range for the tilt becomes 



0.85 < n < 1.16. 
3.2. Cluster Abundance 



(11) 



The next external constraint we consider is the abun- 
dance of clusters of galaxies, which constrains the mat- 
ter power spectrum at interm ediate scales. We adopt 
the empirical fit of Eke et al] (1996) for a flat universe 
er 8 = (0.52 ± 0.08)ft- u oz+u lJS V The value of cr 8 is well 
conver ged w ithin 1 a among different authors (Viana & 
Liddlc 1999; |Pen 1998). This is because the cluster abun- 
dance depends strongly on erg due to its appearance in the 
exponential of a Gaussian in the Press-Shechter formalism. 

We take the amplitude at COBE scales with a 14 % nor- 
malization uncertainty (95 % confidence) together with the 
95 % confidence range coming from the cluster abundance 
to obtain an allowed region that is a function of f2 m , h and 
n and can be roughly described by 



0.27 < n° m 76 hn < 0.35 , 



(12) 



assuming no tensor contribution to COBE and Vl b h 2 = 
0.028. These assumptions lead to the most conservative 
constraints on the (f2 m ,/i) plane. The lower limit comes 
from undershooting erg which is only exacerbated with the 
inclusion of tensors. It also depends on the upper limit on 
n, which is maximized at the highest acceptable baryon 
density £l b h 2 = 0.028. The upper limit comes from over- 
shooting erg and depends on the lower limit on n, which 
only tightens with the inclusion of tensors and lowering of 
the baryon density. 

3.3. Baryon Fraction 

The third external constraint we consider is the baryon 
fraction in rich clusters, derived from X-ray observations. 
The observed baryon fraction shows a slight increase out- 
wards, and the true baryon fraction inferred for the entire 
cluster depends on the extra polation. The est imates range 
from (0.052 ±0.0025)/i" 3 / 2 ( [White fc Fabian| 1995; l owest 
estimate) to (0.076 ±0.008)/t 3 / 2 ( |Arnaud fc Evra7d| 1999; 
highest estimate) for rich clusters. We take the 2cr limits 
to correspond to these two extreme values. We remark 
that very similar constraints are derived from the Sunyaev 
Zcldovich effect for clusters as long as h = 0.5 — 1.0: | My- 



ers et al 



_ (1997) derive (0.061 ± O.Oll)/^ 1 
all (2000) give (0.074± 0.009)/i-\ Adding baryons locked 



and Grcgo et 
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into stars to that in gas inferred by X-ray observations, 
and assuming the cluster baryon fraction represents the 
global value (White et al. 1993), we have fb = Vtb/Q m 



constrained as 

0.052/i~ 3/2 + 0.006/T 1 < f b < 0.076/i~ 3/2 + 0.015/1" 1 

(13) 

This relation is used to convert the constraints on ilbh 2 
into the Q m vs h plane. 

3.4. Age 

We take the lower limit on the age of the universe to be 
to > HGyr based on stellar evolution. While this is not 
based on statistical analysis, no aut hors have ever c laimed 
a cos mic a ge less than thi s value ( Gratton et al. 1997; 
Reid| 1997; Chaboyer et al. 



1998). 



3.5. Allowed Region 

We display all our constraints in the (Q m ,h) plane in 
Fig. 0. 

Combining the range 0.019 < Vt b h 2 < 0.028 from BBN 
and the CMB, together with the constraint on Ob/f2 m from 
the baryon fraction (Oh leads to the range 



0.019 



0.076/i 1 /2 + 0.015/1 



<Q m < 



0.028 



0.052^/2 + O.OO6/1 ' 



(14) 



which is plotted as the solid contours labeled fb in Fig. p0| . 
We also plot the more conservative limits derived from 
taking the 2a extremes of the extreme baryond fraction 
measures (0.076 0.092 and 0.052 0.047) as dashed 
lines. 

We convert the cluster abundance constraint using the 
range in tilts acceptable from the CMB constraints and 
the limit VL b h 2 < 0.028 from BBN (0.85 < n < 1.16) and 
find 

0.15 < n m h 13 < 0.32. (15) 
This range is displayed by the contours labelled "erg" . Fi- 



nally the constraints tu > HOyr and li < 218 are labelled 
as " V' and "li" roapootivoly. 



The shading indicates the parameter space within which 
a model consistent with the CMB and external constraints 
may be constructed. Dark shading indicates the region 
that is also consistent with the stronger nucleosynthesis 
bound of ilbh 2 < 0.023. This does not mean that all mod- 
els in this region are consistent with the CMB data. To 
construct a viable model for a given (fl m ,h) in this region, 
one picks a tilt n between 0.85 — 1.16 consistent with the 
cluster abundance constraint (pL 2|) and then a baryon den- 



sity consistent with H 2 (llCf and n b h 2 < 0.028 (or 0.023) 



In Fig. [j] (solid lines), we verify that the power spectrum 
prediction of a model so constructed is a good fit to the 
data. Here x 2 — 28.5 for the 30 data points to be com- 
pared with x 2 — 27.2 for the model optimized for the CMB 
alone. 

3.6. Consistency Checks 

There are a variety of other cosmological measurements 
that provide alternate paths to constraints in the (n m ,h) 
plane. We do not use these measurements as constraints 
since a proper error analysis requires a detailed consider- 
ation of systematic errors that is beyond the scope of this 




0.2 0.4 



Fig. 10. — Summary of constraints. The shaded region is allowed 
by all constraints considered in this paper. The dark shaded region 
contains models that are also consistent with f2j,/i 2 < 0.023 (sec 
text). 



paper. We instead use them as consistency checks on the 
adiabatic CDM framework. 

Hubble constant. A combined analysis of secondary dis- 
tance indicators gives h — 0.71 ±0.0 4 for an assumed LMC 
distance of 50 kpc ( |Mould ct al] 2000). Allowing for 
a generou s uncertainty in the distance to the LMC (see 
Fukugita 2000 for a review) these values may be multi- 
plied by 0.95 — 1.15 and this should be compared with our 
constraint of 0.6 < h < 0.9. 

Cosmic acceleration. The luminosity distance to distant 
supernovae requires f2 m < 0.48 for flat A models if the 



syste matic errors are no worse than they are claimed ( Hies: 
et al| 1998; |Perlmuttcr et al.| 1999). This limit should be 
compared with our constraint of Q m < 0.6. 

Mass-to-Light Ratio. The J7 m constraint we derived us- 
ing the range 0.019 < fl b h 2 < 0.023 (see Fig. [l(| dark 
shaded region) and the cluster baryon fraction corresponds 
to M/Lb — (350 — 600) h, wh ich is roughly co nsistent with 



M/Lb for rich clusters (e.g. Carlberg et al. 1997). A yet 
larger f2 m (tt m > 0.45) would imply the presence of a sub- 
stantial amount of matter outside clusters and galaxies, 
whereas we h ave some evidence indicating the contrary 
Raiser ct al. 1998). 

Power Spectrum. The shape parameter of the transfer 
function is T ss Vl m h e xp[— 17&(1 — l/f2 m )] « 0.22 — 0.33 
for our allowed region (Sugiyama 1995). This is close to 
the v alue that fits the galax y power spectrum; T = 0.2 — 



0.25 QEfstathiou et al.| 1990; |Peacock fc Doddsj 1994). On 



smaller scales, the Lya forest places constraints on the 
amplitude and slope of the powe r spect rum near k ~ 1 



h Mpc " 1 at z ~ 3 (jCroft et alj 1999; pVlcDonald et a] 



1999). |McDonald et a| (1999) map these constraints onto 
cosmological parameters within ACDM as n = 0.93 ± 0.10 
and cr 8 = 0.68 + 1.16(0.95 - n) ± 0.04. 
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Clu ster Abundance Evolution. The matter density Q 



HI- 



can b ; inferred from evolution of the rich cluster abun- 
dance ( |bukbir fc Blanchard| 1992), but the result depends 
sensitively on the estimates of the cluster masses at high 
redshift. Bahcall & Fan (1998) argue for a low density 
universe fl m = |Blanchard fc Bartlctt| (1998) and 



ulaincn et al. 1999) and lepton asym metry (e.g. Lesgour- 
gues fc Pclosc] 2000; |Esposito ct "at| 2000). On the CMB 



Reichart et al.| (1999) favor a high density Q w 1, while 
Eke et al. (1998) obtain a modestly low density universe 
Vl m = 0.36 ±0.25. 

Peculiar Velocities. The results from peculiar velocity 
flow studies are controversial: they vary from Q m = 0.15 
to 1 depending on scal e, method of analysis and the biasing 
factor (see e.g. Dekel 1999 for a recent review). 

Local Baryons. The CMB experiments require a high 
baryon abundance. The lower limit (together with a mod- 
est red-tilt of the spectrum) is just barely consistent with 
the high baryon abundance option from nucleosynthesis. 
The required baryon abundance is still below the maxi- 
mum esti mate of the bary on budget in the local universe 



0.029/1" 1 QFukugita et al] 1998), but this requires 3/4 of 
the baryons to reside near groups of galaxies as warm and 
cool gas. 

4. FUTURE DIRECTIONS AND IMPLICIT ASSUMPTIONS 

A useful aspect of our approach is that one can ask 
how the allowed parameter space might evolve as the data 
evolves. More specifically, what aspect of the data can 
make the allowed region qualitatively change or vanish al- 
together? If the data are taken at face value, what the- 
oretical assumptions might be modified should that come 
to pass? 

An increase in the precision with which the acoustic 
scale is measured may lead to a new age crisis. It is note- 
worthy that the secondary peaks will eventually provide a 
substantially more precise determination of the scale due 
to sample variance limitations per patch of sky, the multi- 
plicity of peaks, and the effects of driving forces and tilt on 



the fir pt peak [ace Appendix cq. ( A? )]. Indeed, conaiatency 
b etwe bn the determinations of this scale from the various 
peaks will provide a strong consistency check on the under- 
lying framework. If the measurements were to determine 
an equivalent l\ < 200, then to < 10 — llGyrs in a flat 
mology with 9, b h 2 = 0.010; taking n b h 2 = 0.03 dc 



A 

creasd s the age by 1 Gyr and exacerbates the problem 
Such a crisis, should it occur, can only be mildly ame- 
liorated by replacing the cosmological constant with a dy- 
namical "quintessence" field. Because increasing the equa- 
tion of state w from —1 reduces both l\ and the age, only 
a relatively extreme choice of w ^> — 1/3 can help substan- 
tially [see eqn. (0)]. This option would also imply that the 
universe is not accelerating and is in conflict with evidence 
from distant supernovae. However, other solutions may be 
even more unpalatable: a small positive curvature and a 
cosmological constant or a delay in recombination. 

As constraints on the tilt improve by extending the dy- 
namic range of the CMB observations and those on Hi 
by resolving the second peak, one might be faced with a 
baryon crisis. Already Vt b h 2 — 0.019 is only barely allowed 
at the 95% CL. Modifications of big-bang nucleosynthesis 
that allow a higher baryon density for the same deuterium 
abundance are difficult to arrange: current directi ons of 
study include inhomogeneous nucleosynthesis (e.g. Kain- 



side, there are two general alternatives. The first possibil- 
ity is that there is a smooth com ponent that bo osts the 
relative height of the first peak ( Bouchct ct al.| 2000). 
That possibility can be constrained in the same way as 
tilt: by extending the dynamic range, one can distinguish 
between smooth and modulated effects. The direct observ- 
able in the modulation is the ratio of energy densities in 
non-relativistic matter that i s co upled to the CMB versus 
the CMB itself [i?*, see eqn. (A4)], times the gravitational 
potential, all evaluated at last scattering. The second pos- 
sibility is that one of the links in the chain of reasoning 
from the observables to the baryon and matter densities 
today is broken in some way. 

It is noteworthy that there is no aspect of the CMB data 
today that strongly indicates missing energy in the form 
of a cosmological constant or quintessence. An Einstein- 
de Sitter universe with a high baryon density is still vi- 
able unless external constraints are introduced. Under 
the assumption of a flat A cosmology, tight constraints on 
fl m h 3 - 8 from the peak locations and £l m h 2 from the third 
and higher peak heights should allow f2 m and h to be sep- 
arately measured. It will be important to check whether 
the CMB implications for f2 m are consistent with external 
constraints. 

Aside from acceleration measurements from distant su- 
pernovae, the missing energy conclusion finds its strongest 
support from the cluster abundance today through as and 
the cluster baryon fraction. Changes in the interpretation 
of these measurements would affect the viability of the 
Einstein-de Sitter option. 

The interpretation of the cluster abundance is based on 
the assumption of Gaussian initial conditions and the abil- 
ity to link the power spectrum today to that of the CMB 
through the usual transfer functions and growth rates. 
One possibility is that the primordial power spe ctrum has 
str ong deviations from power-law behavior (e.g. Adams ct 
al. 1997). Just like tilt, this possibility can be constrained 
through the higher peaks. 

A more subtle modification would arise if the neutrinos 
had a mass in the eV range. Massive neutri nos ha ve lit- 
tle effect on the CMB itself (podclson et alj 1996; [Ma & 



Bertschinger 1995) but stro ngly suppre ss lar ge scale struc- 



ture through growth rates (Jing et al. 1993; Klypin et al 
1993). A total mass (summed over neutrino species) of 
^2 m Vi = leV would be sufficient to allow an Einstein- 
de Sitter universe in the cluster abundance. One still 
violates the cluster baryon fraction constraint. In fact, 
even for lower fi m one can only find models consistent 
with both the cluster abundance and baryon fraction if 
Y^ m Ui < 4eV. These constraints could be weakened if 
some unknown form of support causes an underestimate 
of the dark mass in clusters through the assumption of hy- 
drostatic equilibrium. They could also be evaded if modi- 
fications in nucleosynthesis weaken the upper limit on the 
baryons. 

5. CONCLUSIONS 

We find that the current status of CMB power spectrum 
measurements and their implications for cosmological pa- 
rameters can be adequately summarized with four num- 
bers: the location of the first peak l\ — 206 ± 6 and the 
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relative heights of the first three peaks Hi = 7.6 ± 1.4, 
H 2 = 0.38 ± 0.04 and H 3 = 0.43 ± 0.07. When translated 
into cosmological parameters, they imply fi m /i 3 ' 8 > 0.079 
(or t < 13 - 14 Gyr), n > 0.85, n b h 2 > 0.019, n m h 2 < 
0.42 for flat ACDM models. Other constraints mainly re- 
flect the implicit (with priors) or explicit use of information 
from other aspects of cosmology. For example, our consid- 
eration of nucleosynthesis, the cluster abundance, the clus- 
ter baryon fraction, and the age of the universe leads to 
an allowed region where 0.6 < h < 0.9, 0.25 < Sl m < 0.45, 
0.85 < n < 0.98, 0.019 < n b h 2 < 0.023. The region is 
narrow, but there clearly are adiabatic CDM models vi- 
able at the 95% CL as exemplified in Fig. |l|. The region 
widens and the quality of the fit improves if one allows 
somewhat higher baryons fl b h 2 < 0.028 as discussed in 
this paper. With this extension the tilt can be larger than 
unity n < 1.16 and £l m as high as 0.6. We note that 
in both cases our limits reflect conservative assumptions 
about tensors and reionization, specifically that they are 
negligible effects in the CMB.0 

The constraints on these and other CMB observables are 
expected to rapidly improve as new data are taken and an- 
alyzed. We have identified sets of observables that should 
provide sharp consistency tests for the assumptions that 
underly their translation into cosmological parameters in 
the adiabatic CDM framework. 

With the arrival of precision data sets, the enterprise 
of measuring cosmological parameters from the CMB has 
entered a new era. Whether the tension between the ob- 
servations that is confining the standard parameters to an 
ever tightening region is indicating convergence to a final 
solution or hinting at discord that will challenge our un- 
derlying assumptions remains to be seen. 
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Steigman, M. Turner, and S. Weinberg for useful discus- 
sions. WH is supported by the Keck Foundation; MF by 
the Raymond and Beverly Sacklcr Fellowship in Prince- 
ton; MZ by the Hubble Fellowship HF-01116-01-98A from 
STScI, operated by AURA, Inc. under NASA contract 
NAS5-26555; MT by NASA grant NAG5-9194 and NSF 
grant AST00-71213. 



where the sound horizon at the last scattering surface s = 
J c s dr] — J c s dt/a with c 2 s = 1/3(1 + R) and R = Sp b /Ap 1 . 
\& is the gravitational potential. The asterisk denotes eval- 
uation at last scattering. Baryons modulate the amplitude 
of the oscillation by shifting the zero point by R^ . The 
result is that the modes that reach maximal compression 
inside potential wells at last scattering are enhanced over 
those that reach maximal rarefaction. Note that this am- 
plitude modulation is not equivalent to saying that the hot 
spots are enhanced over cold spots as the same reasoning 
applies to potential "hills" . 



The oscillator equation (Al) predicts peaks in the angu- 
lar power spectrum at £ m = jiiIa where £a is related to s* 
through its projection on t he sky today via the comoving 
angular diameter distance (Hu & Sugiyama 1995) 



D : 



[i + in(i-n A ) - 086 ] 1+1 - 14(1+,0) 



O EJ-2 ( 1 - n A)- - 76 



2d 



(A2) 



where Q\ refers to the density in dark energy with a fixed 
equation of state w = pa/ PA (w = — 1 for a true cosmo- 
logical constant) and the total density is f2 t = Q m + £1a- 
For convenience, we have defined the dimcnsionless an- 
gular diameter distance d which scales out the effect of 
the expansion rate during matter domination; hence it is 
equal to unity for an Einstein-de Sitter cosmology. More 
specifically, I a = ttD/s* or 



f z*\ 1 ' 2 

172d (w) 



(A3) 



'R* 



In 



y/1 + R,+y/R,+r,R, 



1 + ^fr7RZ 



where the radiation-to-matter and baryon-to-photon ra- 
tios at last scattering are 



r* = p r (z*)/p m (z*) = 0.042w m 1 (z»/10 3 ) , 
R* = 3p b (z*)/4:p y (z*) = SOw^/lO 3 ) -1 



(A4) 



APPENDIX 



with a redshift of last scattering given by 



SCALING RELATIONS 

The phenomenology of the peaks can be understood 
through three fundamental scales which vary with cosmo- 
logical parameters: the acoustic scale £a, the equality scale 
£ eq and the damping scale £jj. 

We begin by employing an idealized picture of the 
photon-baryon fluid before recombination that neglects 
dissipation and time variation of both the sound speed 
c s and the gravitational driving forces. Simple acoustic 
physics then tells us that the effective te mperature pertur- 
bation in the wavemode k oscillates as ( Hu fc Sugiyame] 
1995) 

AT{r],,k) = [AT(0, k) + R^] cos(fcs*) - R*V . (Al) 

2 This assumption is not conservative when considering likelihood 
constraints from the CMB alone. The presence of tensors substan- 
tially weaken the upper limit on n. 



1008(1 + 0.00124^°' 74 )(1 + ciu4 2 ) , 



0.0783w" 



39.5a;"- 76 )- 1 



c 2 = 0.56(1 + 21. W' 8 )- 1 



(A5) 



Here we use the shorthand convention uii, — Q, b h 2 and 
bj m = Q, rn h 2 . Baryon drag works to enhance m =odd over 
m =even peaks in the power. 

These simple relations are modified by driving and dis- 
sipative effects. The driving effect comes from the decay 
of the gravitational potential in the radiation dominated 
epoch which enhances the oscillations and leads to an in- 
crease in power o f approximately a factor of 20 for £ > £ cq 
( Hu fc Sugiyama 1995) where 



4 q = {2Q m H 2 z C(1 ) 1 / 2 D « 438d^ 2 . (A6) 
It also introduces a phase shift to the oscillations such that 



10 



the mth peak of a scale invariant (n = 1 model) is atf] 



= £ A {m - (p) 

r* nO-i 



0.267 



0.3 



(A7) 



Tilt also mildly affects the location of the peaks especially 
the first which is broadened by radiation effects; around 
n = 1 (and fl m h 2 = 0.15), the change is approximately 



Ah 

A£ 2 

£2 
Ah 



0.17(n- 1) , 
Q.033(n- 1) 
0.012(n-l). 



(A8) 



The matter dependence is weak: for Vl m h 2 = 0.25 the 
coefficient 0.17 is reduced to 0.15 for £\. 

The other effect of radiation driving is to reduce the 
baryon drag effect by reducing the depth of the potential 
wells at z* . The baryon drag effect is fractionally of order 
J?**(r7*)/*(77 ini ti a i) w R*T(k) where T{k) is the matter 
transfer function and k is the comoving wavenumber in 
Mpc -1 . The transfer function quantifies the decay of the 
potential in the radiation dominated epoch (see e.g. Eisen- 
stein & Hu 1999 for a fit). The break in the transfer func- 
tion is also given by the horizon scale at matter-radiation 
equality so that it appears on the sky at £ cq . Shifting 
the equality scale to raise £ oq by raising the matter con- 
tent decreases the overall amplitude of the oscillations but 
increases the odd-even modulation leading to somewhat 
counterbalancing effects on the peak heights. 
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Fig. All. — A model power spectrum based on the fundamental 
scales £\, £ e q, 1-q and the baryon-photon ratio R* which modulates 
the amplitude of the oscillations as R,T(£/D) where T(k) is the 
transfer function. 

Finally, the acoustic oscillations are dissipated on small 
scales. The quantitative understanding of the effect re- 
quires numerical calculation but its main features can be 

3 The coefficients are from fits to the first peak at Qf,h 2 = 0.02. 
For better accuracy, replace the coefficient 0.267 with 0.24 for £2 or 
0.35 for £3. Note that the fractional change made by the phase shift 
decreases with m. 



understood through qualitative arguments. Since the os- 
cillations dissipate by the random walk of the photons in 
the baryons, the characteristic scale for the exponential 
damping of the amplitude is the geometric mean between 
the mean free path Ac = (x e 7i e 0T<x) and the horizon 
scale 



-l/2 z l/2 



[vT 



(A9) 



— 1/2 

under the Saha approximation x e cx ui b so that kp 



(77* Ac) 



-1/2 



5/4 1/4 1/4 



CX 2* Ul 



Numerically, the scaling is 



slightly modified to (refitting values from Hu & White 
1997) 



ID 



k D D: 



2240d 



5/4 



,0.24, 



-0.11 



[(1+n) i/2_ r i/2 ]1/2 V10007 m 

(A10) 

Compared with the acoustic scale I a, it has a much 
stronger dependence on and the redshift of recombi- 
nation z» . We show a model spectrum^ obtained this way 
for a ACDM cosmology with fl m = 0.35, O t = 1, h = 0.65, 
uib = 0.02, and w = — 1 in Fig. All] . For this cosmology, 



the three fundamental scales are £ cq = 149 (£\ = 221), 
£a = 301, and £n = 1332. The dependence of the mor- 
phology of the acoustic peaks on cosmological parameters 
is controlled by these three scales. Around the fiducial 
ACDM model with the parameters given above 



At, 



-O.llAi 



0.24- 



0.07 



Auj b 



-0.17—^-1.1- 



fit 



Ah Auj b 
-O.llAto - 0.48— + 0.07 

h u>b 



0.15- 



-1.4- 



.Afit 
fit 



(A15) 



where the leading order dependence is on f2 t and h 



AL 



O.llAw + 0.5- 



0.17- 



AO/ 



fiA 



1.1- 



Afit 
fit 



-O.llAw 



Ah 
~h 



0.59- 



A^„ 



1.4^, (Ale, 



rum is obtained by following a construction 
Whitel (1997): the damping envelope is 

V e =e X p[-(£/£ D ) 1/2 ], (All) 

yielding acoustic oscillations of the form 

V t Ai = [1 + R.*T(£/D)]V e cos[n(£/£ A + 0)] - R*T(l/D) ; (A12) 

the potential driving envelope is 

Vl fa l + 19exp(-1.4< eq /*). (A13) 

The spectrum is then constructed as 



(AT,) 2 cx 



V t V\ 



1 + (£ A /2£)S 



(A14) 



where we have added an offset to the oscillations to roughly account 
for projection smoothing and the Doppler effect and forced the form 
to return to Vt above the first, peak to ac count for the early ISW 
effect from the radiation (Hu &; Sugiyama 1995). This mock spec- 
trum should only be used to understand tlie qualitative behavior of 
the spectrum. 
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which depends more strongly on Q m , and 



A£ D 



-O.llAtu - 0.21 



0.20 



: D 



Ato b 



-0.17 



Afl A 



1.1- 



AOt 



Ah 

-O.llAw - 0.42 

h 



0.20 



Auj b 



-0.12- 



AO- 



Or 



1.4- 



(A17) 



which depends more strongly on the baryon abundance 
Note that the sensitivity to fit incre ases from the often 
quoted — Q.5Af2t/fit as J7a increases ( Weinberg 2000; M. 
Turner, private communication). 

Ideally one would like to extract these three numbers 
and the baryon-photon ratio i?* directly from the data. 
The acoustic scale is readily extracted via the position of 
the first and/or other higher peaks. The other quantities 
however are less directly related to the observables. We 
instead choose to translate the parameter dependence into 
the space of the observations: in particular the height of 
the first three peaks. 

The height of the first peak 



H 



( AT lx 
VATio 



(A18) 



may be raised by increasing the radiation driving force 



(lowering 



eq 



or the baryon drag (raising w&j. How- 



ever it may also be lowered by filling in the anisotropics 
at £ w 10 through the ISW effect (raising Qa or w, or 
lowering fit), reionization (raising the optical depth r), or 
inclusion of tensors. Each of the latter effects leaves the 
morphology of the peaks essentially unchanged. Because 
Hi depends on many effects, there is no simple fitting for- 
mula that describes it. Around the ACDM model with 
Hi = 7.4, it is crudely 



AHi 



Alu,, 
-0.5 

2.5An 



0.4^^ - 0.5Afi A + 0.7Afi t 



0.02 

IAt - 0.3Aw - 0.76 



Ar 



1 + 0.76r 



where the tensor contribution r = 1A(AT^ /AT$) 2 . 
This scaling should only be used for qualitative purposes. 

The height of the second peak relative to the first is 
written 



dependence on the matter comes from two competing ef- 
fects which nearly cancel around the ACDM: increasing 
uj m (lowering £ eq ) decreases the radiation driving and in- 
creases H2 but also increases the depth of potential wells 
and hence the modulation that lowers i?2- 

For the third peak, these effects add rather than cancel. 
When scaled to the height of the first peak, which is also 
increased by raising the baryon density, the flbh 2 depen- 
dence weakens leaving a strong dependence on the matter 
density 



Ha 



( AT t . 



\AT ei 



2.17 



1 



{—)' 
V0.044/ 



0.044/ 
1 + 1.63 ( 1 



^°„ 59 (3.6)" 

-1 



0.071 



(A21) 



where £3/^1 ~ 3.6. 
where H 3 = 0.50 



Around the fiducial ACDM model 



AH 3 



0.4lAn- 0.31 



Aiv b 

UJ b 



0.53- 



(A22) 



We emphasize that the phenomenology in terms of £a, 
£ eq and £d is relatively robust, predictive of morphology 
beyond the first three peaks and readily generalizable to 
models outside the adiabatic cold dark matter paradigm. 
The specific scalings of Hi , H2 and H3 with cosmological 
parameters are only valid within the family of adiabatic 
CDM models. Furthermore, as the data continue to im- 
prove, the fits must also be improved from their current 
few percent level accuracy. The number of phenomenolog- 
ical parameters must also increase to include at least both 
the heights of the peaks and the depths of the troughs for 
all observed peaks. 



Hi 



[An 



\AT ei 



0.925w^ 18 (2.4) r 



[1 + (wb/0.0164) 



12wO l 52 ll/5 



(A19) 



where n is the scalar tilt and £2/^1 ~ 2.4. This approxi- 
mation breaks down at high iOb and uj m as the second peak 
disappears altogether. In the ACDM model with n = 1, 
H2 = 0.51 and parameter variations yield 



AH 2 
H 2 



0.88An-0.64 



A^ 

UJb 



0.14- 



(A20) 



The effect of tilt is obvious. Baryons lower H2 by in- 
creasing the modulation that raises all odd peaks. The 
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